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Abstract 

Surfaces which do not exhibit surface melting below the melting point (non- 
melting surfaces) have been recently observed to sustain a very large amount 
of overheating. We present a theory which identifies a maximum overheat- 
ing temperature, and relates it to other thermodynamical properties of the 
surface, in particular to geometrical properties more readily accessible to ex- 
periment. These are the angle of partial wetting, and the nonmelting-induced 
faceting angle. We also present molecular dynamics simulations of a liquid 
droplet deposited on Al(lll), showing lack of spreading and partial wetting 

in good agreement with the theory. 
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For a long time crystal overheating above the bulk melting temperature Tm has been 
believed to be impossible, at least in the presence of a free clean surface. The standard 
argument is that surface premelting will always take place and act as an ubiquitous 
seed for the liquid to grow. The well-known surface melting of Pb(llO) provided a 
first microscopic evidence of how liquid nucleation takes place on a solid below T^. It was 
only a little later that simulations of Au(lll) |^ and newer experiments on Pb(lll) and 
Al(lll) [0 demonstrated microscopically that the opposite could also happen, namely that 
certain surfaces may exhibit nonmelting up to and in fact even above the melting point [^H. 
A solid bounded by such surfaces can therefore be overheated, although in a metastable 
state, above T^. Metois et al. have first shown that small Pb particles with strictly (111) 
facets are easily overheated by a few degrees above 0. Even more strikingly, Herman 
et al. found that a flat nonmelting Pb surface can be overheated by as much as 120 K 



above Tm [|T0|- This implies that the free energy of a crystal surface can have a local 
minimum for zero liquid thickness. As in other nucleation problems one should thus expect 
the metastable overheated state to survive up to some instability temperature Tj > T^, 
where the barrier finally disappears (fig. |I], inset). At present, however, there is no further 
available understanding of this phenomenon. In particular, there are no means to calculate 
Tj and possibly connect it with other quantities which are more readily measurable in a 
surface experiment. At a more microscopic level, it is very desirable to understand the 
different behavior of a nonmelting and of a melting surface, against nucleation of the liquid. 

In this Letter, we introduce a simple theory of surface nonmelting which predicts the 
existence of a Tj, and connects its value with apparently unrelated geometrical quantities. 
These are the partial wetting angle 6m which a drop of melt will form with that crystal 
surface at T = T^, and the faceting angle 6c of a vicinal surface. The angle 6^ has also been 



rather commonly measured in the past, a few early examples being the (0001) face of Cd |Tl 



and the (100) faces of several alkali halides The nonmelting-induced faceting |T3||n 
angle 6^ has been well characterized experimentally and theoretically for (111) vicinals of 
Au |[13|JT5[| , Cu [|l6l and Pb |]l3|,|l4|,|l^|l8l . The connection we find between 6m, 6c and 



Ti offers new insight in nonmelting surfaces. At a microscopic level, we substantiate this 
connection with molecular dynamics (MD) simulations of Al(lll), which demonstrate both 
the non-spreading of a liquid drop at Tm, and the overheating of the fiat face. The predicted 
relationship between 6m, dc and Tj is found to be in excellent agreement with the simulation 
results, as well as with experiments. 

(i) Theory: Consider a liquid film of thickness sandwiched between semiinfinite solid 
and vapor, and let d. grow from zero (no liquid) to a finite value. The change in free energy 
per unit area takes the standard form 

AF(£) = pU{l - T/Tm) + A7(£) (1) 

where p is the liquid density, L the latent heat of melting, and A7(£) the difference between 
the overall free energy of the two interacting solid-liquid (SL) and liquid- vapor (LV) interfaces 
separated by a distance and the free energy of the solid- vapor (SV) interface. By definition, 
A7(0) = 0. Assuming short-range forces only, this term can phenomenologically be written 
as A7(^) = A7oo[l — exp(— where A700 = 7sl + 7lv — 7sv is the net free energy 
change upon conversion of the SV interface in two non-interacting SL and LV interfaces, 
and ^ is a correlation length in the liquid. For a melting surface A700 < 0, and, for 
Tu, < T < Tm, AF will have a minimum at £o(^) = ^ In [T^l A7oo|/(rm — r)Lp^] which 
is the mean-field thickness of the melted film The wetting temperature defined by 
4(T^) = is = r„ (1 - |A7oo|/LpO- 

For a nonmelting surface, A700 > 0, and we move over to T > Tm- Here, AF(£) will 
instead have a local minimum at £ = 0, the absolute minimum for £ ^ 00, and a maximum 
at a critical thickness 



4(T)=ein 



^mA7o 



(2) 



{T-Tm)Lpi 

as shown in fig. |1|. The local minimum at £ = signifies metastability of the crystalline 
surface for T < Tj, the maximum overheating temperature. The minimum disappears when 
= 0, yielding 



Above Ti, the crystal surface will melt, no matter what its initial state is. In particular 
a surface which is initially crystalline will wet itself with a liquid film, which will grow, 
and gradually melt the whole crystal. Hence Tj can be seen as a non-equilibrium wetting 
temperature, or, more accurately, as a spinodal point for the overheated solid surface. For 
Tm < T < Ti, the predicted behavior is that typical of a nucleation problem. If the surface is 
prepared initially with a melted film of thickness i > ic (upper vertical arrow in fig. [I|), then 
melting will proceed, and i will grow to infinity, reaching full equilibrium. If, conversely, 
the starting thickness is less than i^, then the surface will recrystallize, to reach the local, 
metastable minimum at £ = (lower arrow in fig. This peculiar behavior was first found 
and described in detail in an early simulation of the nonmelting surface Au(lll) P|. 

We now show that there is a simple connection between Tj, and the macroscopic non- 



wetting angle 6*^ at T = Tm- Following Nozieres [ffH], the angles 6'lv, ^'sl, formed by a drop 



of melt onto a nonwetting surface of the same material (fig. 0), satisfy the equations 

7sv = 7lv cos 6'lv + 7sl cos 9sl (4) 

i?Lv sin 6'lv = ^SL sin 6'sl (5) 

where -Rlv, -Rsl are the radii of respectively the LV and SL drop boundaries (suppos- 
edly spherical). Eq. (^) is simply the balance of lateral forces, while eq. (^ follows 
from simple geometry. Laplace's pressure equation P = 2'y/R determines the shape ra- 
tio x(T) = Ru//Rsh = sin6'sL/ sin6'Lv = [7LV-PsL(r)]/[7sL-PLv(r)]. Since Psl oc (T - T^) 
near Tm, we expect 6'sl to switch from negative for T < to positive for T > Tm- At 
T = Tm, X = 6'sl = 0, i?LV = oo, and 6'lv = 6'm where 

cos6'm = 1 (6) 

7lv 

Comparison of (^ with (|^) shows that knowledge of A700 at T = Tm determines both Ti and 
dm, which are monotonically related by 



Ti 27lv . 2 

—— = 1 + ^ ^ sm — . (7) 
For a nonmelting surface there is a second important angle Oc-, which is the nonmelting- 



induced faceting angle |[T3|Jl^ . Consider vicinal faces tilted at an angle 9 away from the 
nonmelting face. At T = Tm there are two well-defined free energy minima (solid, £ = 
and liquid, £ = oo). We can thus draw [|1^ the two projected surface free energy branches 
a{9) = '~f{6)/ cos6 as a function of the step density t = |tan0|. The two branches are 
approximately given by the standard expressions 

(Ts{e) = -fsv + l^t + gt^ (8) 



{0) = {1ly + 1sl)ViT¥ (9) 



where /i and g are the step free energy and the step-step repulsion on the solid surface. 
Here we have further assumed that 7sl is approximately independent of 6. The faceting 
angle is given by 9c = arctantc which satisfies the double tangent construction: cq -|- cito = 



Isv+fJ'to+gtl, ci = ix+3gtl, and Co+Citc = (7lv + 7sl) y^l + t^, ci = (7lv + 7sl) tc/\J^ + tl- 
A particularly simple solution is obtained if the cubic (step-step repulsion) term gt^ can be 
ignored, whence to = 0, cq = 7sv, ci = [(7lv + 7sl)^ ~7sv]"'^^^ (note the nonanalyticity of as 
at t = 0) and tc = {[(7lv + 7sl) /7sv]^ — l}^^"^. Even when this approximation cannot be 
made, and to is nonzero (as is the case for Pb(lll) where arctanto — 2°), the above is 



still a pretty good approximation to the faceting angle, which is therefore simply related to 
A700: 

cos^^,^fl + ^^ . (10) 
V 7sv / 

Eq. (H) shows that 9^ is identical to both the droplet angles 6'lv, 9sl at a single temperature 
Tu > Tm, satisfying x(T„) = i?Lv/-RsL = 1- We note that 9c is slightly smaller than 9m- The 
outer droplet angle 6'lv will therefore decrease from 9m to 9c to zero when T is raised from 
Tm to Tu to Ti. Finally, we observe that a physical upper bound for 9m and 9c is given by 
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A7oo ^ 7lv, whence 9m ^ 90°, and 9c ^ 60°, i.e., a melt must at lest partially wet its own 
solid. 

(ii) MD simulations: Choosing Al(lll) as our test case, we have simulated its behavior 
at and above Tm using the recent accurate glue potential of Ercolessi and Adams |^ , derived 



by fitting to first-principles data. First, the approximate bulk melting point for this potential 



was determined using the phase coexistence technique |^ and found to be Tm = 939 ± 5 K 
(against an experimental value of = 933.6 K). Then a 16-layers slab with 3 rigid 
bottom layers, one free surface, 224 atoms per layer and x-y periodic boundary conditions, 
was studied as a function of T. As Tm was reached and crossed, the surface remained 
crystalline (metastable) as expected, up to a large Tj = 1088 ± 18K = Tm + (149 ± 18) K, 
even for very long (2 ns) runs. On the basis of our theory, using the known values of 
p = 0.0534 A'^ L = 105.4 meV/atom [gg, and an estimated ^ = 2.6 ± 0.3 A [|2 



we 



extract from (^ A700 = 2.3 ± 0.4meVA~^. Inserting in eq. (^, with an estimated value 
of 7lv = 46.6 meVA^^ (obtained with a separate simulation of the free liquid surface at 
T = Tm) we finally predict 9m = (18 ± 2)° and, with a value 7sv = 54.3meVA~^ [|20| , 
= (16±2)°. 

To check this prediction, we have prepared an 861 Al-atom cluster which is fully melted 
and forms a liquid drop already at 900 K By depositing this Al drop on any given 

Al surface, we can learn about its wetting habit. We deposit it first on the Al(llO) face, 
which is prone to melting 0. At T = 930 K (below Tm, but above T^), the drop spreads out 
completely within 100 ps (fig. |^a-c). However, when deposited on the nonmelting Al(lll) 
face it does not diffuse away, but rather settles down as expected with well-defined exterior 
and interior angles whose azimuthal average (9) we can extract. By increasing temperature 
across T„, from 930 to 945 K, we find that (^lv) changes from (24 ± 3)° to (21 ± 1)°, and 
(6'sv) from — (6'lv) (the droplet is essentially crystallized) to (44 ± 6)°. By interpolation we 
extract 9m = (22 ± 3)°, in fairly good agreement with the predicted value (18 ± 2)°. The 
approximate values of ^, and of 7lv (about 20% lower than its experimental value, with 
our potential) constitute sources of error. Additional discrepancies are to be attributed to 



the macroscopic and phenomenological nature of the theory, which should in principle be 
improved to include fluctuations and flnite-size effects. On the simulation side, one could 
consider in the future flnite-size scaling as a possibility. 

(Hi) Connection with experiments: We are not aware of measurements of 6m, dc or Ti 
on Al(lll), for which we have thus a direct prediction. On Pb(lll) van Pinxteren and 
Frenken measured 6c = (14.7 ± 1.4)°, whence using eq. (0), and for 7sv — 34meVA~^ 
we obtain A700 = 1.2 ± 0.2 meVA~^. Furthermore, using p = 0.033 A~^, L = 50meV/atom, 
and ^ ~ 2.7 A (averaging data from ref. as suggested in fl^), we obtain via eq. (H) 



Ti = Tm+(150±30) K. This is in rather good agreement with the experimental result, T^^^ ~ 
T^P + 120 K [0]. From eq. (|), using 7lv ^ 28meVA-2, we also predict 9^ = (16 ± 1)° 
for a droplet on Pb(lll) at T = T^. For Al(lOO), another nonmelting surface p5[, we flnd 
by simulation Tj = 1025 ± 5 K, and assuming ^ ~ 3A we predict 6^ — 15° and 6^ ~ 13°. 
For Cd(OOOl), where 6^ = (37 ± 1)° and using p = 0.043 A'^, L = 64meV/atom, 
7lv = 40 meVA~^, and again a guessed ,^ ~ 3 A, we get Tj ~ + 580 K, close to twice the 
melting temperature (594 K). Application of this scheme to Ge(lll) or NaCl(lll) appears 
instead problematic, due to the essential role of long-range forces in these cases. In fact 
Ge(lll) has a negative Hamaker constant, which is probably related to its nonmelting 
behavior |2^, while long-range Coulomb forces are likely to be relevant to the nonmelting 
of NaCl(lOO). Finally, the present scheme is probably also inapplicable in its simplest form 
to surfaces such as Pb(lOO) or Au(lOO) which undergo incomplete melting |[T7|j27| -p9| . 

In summary, we have described new results on the nonmelting crystal surfaces. A simple 
theory of nonmelting is given, which describes the metastable solid surface above T^, up to 
a maximum overheating temperature Tj which acts as a spinodal point. This temperature 
is found to have a simple connection with the partial wetting angle of the surface by its 
own melt at T ^ T^, and with the tilting angle of melted regions in vicinals undergoing 
nonmelting-induced faceting. Computer simulations on Al(lll), as well as available data on 
Pb(lll) are in good agreement with this theory. 
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FIGURES 

FIG. 1. Critical liquid thickness of a nonmelting surface vs. temperature above (schematic). 
A system with a liquid film thinner than the critical value, will recrystallize for any T between Tm 
and Tj. One with a thicker film will melt completely. Inset: free energy change upon conversion of 
a film of thickness I from solid to liquid. Prom to Tj the solid surface is a local minimum. 

FIG. 2. Shape of a drop of melt onto a nonmelting solid surface of the same material. The two 
interfaces separating solid and liquid (SL), and liquid and vapor (LV) are assumed to be spherical, 
with radii i?sL and ii!LV aiid contact angles ^SL and ^lv respectively. 

FIG. 3. Evolution of an 861-particles liquid drop of Al on substrate of the same material. 
Left column: drop on a surface undergoing surface melting (Al(llO) at T = 0.99 T^). (a) before 
contact; (b) after contact, the drop spreads readily; (c) the drop has been fully absorbed. Right 
column: drop on a nonmelting overheated surface (Al(lll) at T = 1.01 T^). (d) before contact; 
(e) after contact: the drop settles, but does not spread; (f) final drop shape. Darkness of atoms is 
proportional to their square displacement in the run. 
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